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For the non-perturbative analysis of the topological content of Yang-Mills theories, it is essential 
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proposed to use the adjoint modes of the Dirac operator for this purpose. In this talk we analyse an 
implementation of this idea that associates two Weyl fermionic modes in the adjoint representation 
to every gauge field configuration. The densities of these modes provide a filtered image of the 
self-dual and anti self-dual parts of the gauge action density. We present successful tests on the 
performance of this proposal on a set of initial gauge field configurations. 
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1. Introduction 

In the past years several authors have devoted some effort to the analysis of the topological 
structures present in pure Yang-Mills theory on the lattice. Apart from global information such as 
the topological susceptibility, it is interesting to extract local one, such as the instanton size distri- 
bution. In some cases one might be able to check certain ideas and proposals relating these topo- 
logical structures to chiral symmetry breaking and confinement. In this program the battlehorse 
is the roughness of lattice Monte-Carlo configurations, having its origin in divergent ultraviolet 
fluctuations. To tame this noisy background several ideas have been proposed. In particular, cool- 
ing and smearing algorithms produce smoother configurations which are assumed to preserve the 
long-range structure of the original one. However, these methods have been criticised since they 
may produce some distortions on the shape distribution of local structures. A proposed alternative 
is to use filtering methods based on the Dirac operator or other differential operators ^ (see for 
example [Q] and references therein). For the case of the Dirac operator, the main idea is the rela- 
tion between fermions and topology given by the Atiyah-Singer index theorem and the correlation 
between the gauge action density and the local density of the eigenstates of the Dirac operator. In 
practise, however, these methods are not totally free of ambiguities. For example, when fermions 
in the fundamental representation are used to reconstruct the topological charge density the ability 
to reproduce the topological structure depends, in a rather strong way, on the number of modes in- 
cluded in the reconstruction. In any case, even for noiseless configurations the shape of the filtered 
density does not coincide with the corresponding action density. 

In this talk we analyse an alternative proposal, presented in [0], based on the use of the adjoint 
representation of the Dirac operator. The advantage of this method is that it is based upon a single 
mode and gives a perfect match for classical solutions of the equations of motion. In the next 
section we will give details of the method and its lattice implementation. In Section 3 we will test 
its performance by analysing its filtering capacity for a series of initial gauge field configuration. 

2. Description of the method 

For every gauge field configuration we will construct two associated classical Weyl fermionic 
fields in the adjoint representation, \\f±(x), whose densities, ||y/±(x)|| 2 , correspond to a filtered 
image of the self-dual and anti self-dual parts of the initial gauge field action density. The Weyl 
spinorial fields have components Ya±(x) where a is a 2-component spinor index, and a is the 
colour index taking N 2 — I values. From now on we will refer to the fields i//± as adjoint filtering 
modes or AFM. 

An important property that any filtering method should satisfy is that for smooth configu- 
rations the procedure must reproduce the classical structures without distortion. The main idea 
behind the method introduced in Ref. [J|] is that for configurations that are solutions of the classical 
equations of motion (classical solutions) there is an optimal choice of the spinorial field in the ad- 
joint representation that reproduces exactly the shape of the (anti)self-dual part of the gauge action 
density. This is obtained choosing the AFM, \j/±, as the chiral components of the supersymmetric 
zero-mode of the Dirac operator, defined by: 

¥=^Fnv[Yti,7v]V, (2.1) 
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where V is an arbitrary constant spinor. This mode has three interesting properties: 

1. It is a zero-mode: ifoxff = 0. 

2. It satisfies the reality condition, lm(\irf ± (x)) = Va,jc. 

3. Its density is exactly proportional to the (anti)self-dual part of the gauge action density. 

Thus, for a classical solution, if we select a mode satisfying properties 1 and 2, then by virtue of 
the third property, we obtain an optimal image of the field density structure. 

For a general gauge field configuration, our method consists in finding the Weyl-spinor fields 
in the adjoint representation that satisfy properties 1 and 2 as much as possible. Different versions 
of the method follow from the relative importance given to both properties. Previous tests @ Q] 
were done by looking first at the subspace of low-lying eigenmodes of the adjoint Dirac equation 
and then selecting the combination which best satisfies the reality condition as the AFM mode. 

Here we will use an alternative more elegant proposal, also presented in Ref. [||], in which the 
reality condition (property 2) is imposed exactly. The other condition is implemented by requiring 
that the Weyl spinor is an eigenvector of lowest eigenvalue of —0 . For a left-handed (positive chi- 
rality) spinor this becomes —DD, where D = D^o^ (D = D^o^) are the Weyl operators. Actually, 
in the adjoint representation all eigenvalues are doubly degenerate due to euclidean CP invariance. 
Then, using any eigenvector y and its charge conjugate \j/c = 02 V* one can form a colour vector 
of quaternionic matrices 

( ¥l¥c ) = ^ ± (x)a tl , (2.2) 

where = (I, —if), cfy = (I,/?) and T, are the Pauli matrices. The eigenvalue equation becomes 
an equation acting on quaternionic matrix fields, and the reality condition amounts to *P± (x) = 0. 
In summary, the positive chirality AFM mode is defined by the eigenvalue condition: 

-D&V^xjOi = X^ + {x)Gi , (2.3) 

or equivalently 

Op 3 + {x)=W i + {x), where 0+ = -D^Dyrj^fjf , (2.4) 
and 17 and t) are the 't Hooft symbols. The last expression can be expanded and we get 

0± = -8iPl - etj^D^Dy = -SijD* + ie ijk (E k + B k ) . (2.5) 

Notice, that the + operator is a positive definite, real, symmetric operator, so that its eigenvalues 
are positive real numbers. For gauge configurations which are solutions of the classical equations 
of motion the lowest eigenvalue vanishes and generically is non-degenerate. The corresponding 
eigenvector is the supersymmetric zero mode °< (Ej+Bi). 

For arbitrary gauge configurations the minimum eigenvalue of + is non-zero. Its corre- 
sponding eigenvector is, by definition, the AFM mode. Its density provides the filtered version of 
the self-dual part of the action density. The usefulness of the method depends on its capacity to 
eliminate high frequency noise without altering the shape of smooth structures. In this talk we will 
present the results of our tests done on both smooth and rough configurations. 

The previous formulae can be repeated for the negative chirality mode, which provides a fil- 
tered version of the anti-self-dual part of the action density. The corresponding operator is now 
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replaced by 




dijD 2 -ie ijk (E k -B k ). 



(2.6) 



In order to obtain a lattice implementation of the filtering procedure it is more convenient to 
work with the overlap Dirac operator and construct the hermitian, positive definite matrices H\ = 
P±(y 5 D ov ) 2 P±, where P± are the projectors onto the positive and negative chirality modes. The 
two-fold degeneracy associated with CP invariance also holds on the lattice, so that the operators 
can be taken to act on quaternionic vectors. Imposing the reality conditions leads to the 



where the action of Pq on a quaternionic field is given by Pq^ 1 = *P — ^GoTr(ao x i J ). 

The AFM modes are the eigenvectors of lowest eigenvalue of Of and can be obtained by a 
conjugate gradient algorithm. Details of the technique used to compute the overlap Dirac operator 



3. Testing the filtering method 

In this section, we present the results of our tests of the filtering procedure when applied to 
various lattice configurations. First, we apply it to an instanton configuration. Being a classical 
solution the method should work well, but it will allow us to quantify the finite volume and dis- 
cretisation effects. Next we apply it to a series of instanton-anti-instanton (IA) pairs with varying 
separations. This is intended to monitor possible distortions and problems which could occur when 
applied to smooth configurations which are not classical solutions of the equations of motion. Fi- 
nally, we will go back to the single instanton case and add stochastic noise to it. The filtering 
method should be able to reduce this noise and produce a neat image of the underlying smooth 
configuration. 

3.1 Results for Classical solutions 

The first test is done over a set of smooth, Q = 1, SU(2) instanton configurations generated 
by cooling. In this case, the + operator should have one zero-mode corresponding exactly to the 
supersymmetric zero-mode. As will be discussed below, discretisation and finite volume effects 
may shift the corresponding eigenvalue from zero. Still, we observe that the density of the lowest 
eigenvector of the + operator reproduces to an excellent degree the instanton action density. 
Fitting both the action and the AFM mode densities to the continuum formula, we extract instanton 
positions and sizes that differ at most in AX = 0.01 a, and in Ap = 0.05 p, respectively. 

In contrast to the situation in the continuum, the lowest eigenvalue of is different from 
zero. This is due to discretisation errors and finite volume effects. The latter arise for the case of 
periodic boundary conditions (PBC), because Q = l classical solutions do not exist on a periodic 
torus This is, however, not the case if twisted boundary conditions (TBC) are used instead. 

To explore both effects we have generated a large set of SU(2) instanton configurations with 
varying sizes for both PBC and TBC. Fig. [Ta| displays the lowest eigenvalue (X\) of the + operator 
as a function of the inverse instanton size squared a /p . As expected, X\ approaches zero in the 
continuum limit as a 2 /p 4 , with a coefficient of approximately 8 • 10~ 4 . 



Of=P HlP 



(2.7) 



and other numerical aspects can be seen in Ref . [||] 
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Figure 1: (a) The lowest eigenvalue of the + operator, on a size p instanton background, is 
displayed as a function of a 2 /p 2 for two different lattice sizes and boundary conditions. The finite- 
volume effects become sizable for periodic BC and L < 4p. They are negligible for time-twisted 
BC (k= (1, 1, 1), m = (0,0,0)). (b) The three lowest eigenvalues of the + operator, on an IA pair 
background, are displayed in terms of the IA distance (d). Level crossing appears at d ~ 1.8 p. 



In what concerns finite-volume effects, a clear distinction can be observed between periodic 
and twisted boundary conditions, as expected. In the twisted case, classical solutions exist for all 
torus sizes and we observe no corrections to the lowest eigenvalue other that those associated to 
discretisation effects. For periodic BC, however, we observe a large deviation from the expected 
behaviour for L < 4p. 

3.2 Smooth non-classical configurations 

According to our proposal, even for non-classical configurations the lowest eigenvector of 
the + (0~) operator should provide the filtered action densities of the self-dual (anti-self-dual) 
part of the gauge field. Since only for classical solutions there is an exact supersymmetric zero- 
mode, it is not guaranteed that the filtered image of non-classical configurations has no distortions. 
In order to test this, we have analysed a set of instanton-anti-instanton (IA) configurations with 
varying separation. The set is generated by cooling an initial well-separated pair and monitoring 
the different stages of the IA annihilation process. Fig. || displays one snapshot corresponding to 
IA distance d = 3.5 p. We compare the AFM densities with the self-dual and anti-self-dual parts 
of the action density. The agreement is excellent and, as expected, the O lowest eigenvectors 
are only sensitive to objects with the appropriate chirality. In order to obtain a more quantitative 
comparison, we extract the size parameter of the (anti-) instanton by fitting both the AFM and 
action densities. Both determinations differ by 5% at most as long as the IA separation is larger 
than 2p. 

Although the Dirac operator does not have zero-modes on IA backgrounds, we expect the 
lowest eigenvalue in each chiral sector to go to zero as the I-A separation is increased. To show that 



this is indeed the case, we display in Fig. lb the three lowest eigenvalues of the + operator versus 
the IA separation. The smallest one is well described by an exponentially decreasing function of the 
IA separation as d 2 X\ = 16 exp(— 2d/p). For distances d ~ 2p we observe a level crossing in the 



spectrum of the O operators. For even smaller separations the lowest mode no longer reproduces 
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Figure 2: For an instanton-anti-instanton pair, at separation <i = 3.5p, we display 2-d slices of: 
(Left) the self-dual (red) and the anti-self-dual (blue) parts of the action density. (Right) the AFM 
density corresponding to the positive (red) and the negative (blue) chiralities. 




Figure 3: Comparison of the gauge action density (left) and the AFM density (right) for a configu- 
ration generated by applying ten heat-bath sweeps (j3 = 4) to a smooth instanton configuration. 



the IA gauge action density. Hence, we conclude that the operators can identify the components 
of the IA pair as long as d > 2p . 

3.3 Configurations with stochastic noise 

The most important test, prior to its application to Monte Carlo configurations, is to show that 
the procedure does indeed filter out high frequency noise from the starting configuration, exhibit- 
ing its long-range structures. For this purpose, we began by generating several smooth instanton 
configurations and added random noise to them. The way this was done was by applying to the 
configuration a small number of heat bath updates corresponding to a Wilson action with different 
values of /3. The small number of updates and the large values of /3 guarantee that the instanton 
is not destroyed in the process, but considerable noise is added. The results presented here corre- 
spond to an initial instanton of size p = 3.4a, on a periodic/twisted lattice of size 14 4 , to which ten 
heat bath sweeps (with j8 = 30, 20, 8, 7, 6, 5, and 4) have been applied. The process is repeated with 
ten different initial random seeds giving rise to a heated instanton ensemble for each j8 value and 
BC. A characteristic example is shown in Fig. ||, where we display the AFM and the gauge action 
densities for one rough configuration with j8 = 4 and PBC. While the action density dramatically 
roughens under heating, the AFM density is practically insensitive to the high frequency modes. 
The same is observed in the twisted case and for all the j8 values analysed. This is a good proof of 
the extraordinary filtering capacity of our method 

The ability of the method to recover the initial instanton structure is related to the fact that, 
despite the eigenvalue becoming non-zero after the addition of noise, there seems to be no level 
crossing and the AFM is always cleanly separated. This is clearly seen in Fig. || where we display 
the four lowest eigenvalues of the + operator as a function of /3" 1 , for periodic (left), and twisted 
(right) boundary conditions. Notice that not only is there no level crossing, but also the gap seems 
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Figure 4: We display the four lowest eigenvalues of the + operator for the ensemble of heated 
instanton configurations described in section 3.3. For each value of j8 we plot the eigenvalues 
averaged over the ensemble of heated configurations, with error bars corresponding to the variance 
of the sample. The continuum lines correspond to a second order degree polynomial fit in jS 1 . 



to remain constant when the size of the noise (controlled by /3 x ) increases. The data for the lowest 
(AFM) and first excited eigenvalues of O^ 1 are well fitted by a second degree polynomial in j8 _1 , 
which is the form predicted by perturbation theory around the instanton configuration. The constant 
term is fixed to the value obtained for the instanton before heating (which for the ground state is 
determined by finite-volume and finite lattice spacing errors). The remaining coefficients turn to 
be very similar for both eigenvalues and boundary conditions. 

4. Conclusions 

In this talk we have analysed the results of the filtering method proposed in Ref. ^ when 
applied to different lattice configurations, some smooth and some rough. The results are quite 
encouraging and point out to the main difficulties one might encounter when applying the method 
to thermalised Monte Carlo configurations. Additional details and tests can be found in Ref. [|5[ [7]]. 
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